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We develop the general scheme for modified f(R) gravity reconstruction from any realistic FRW 
cosmology. We formulate several versions of modified gravity compatible with Solar System tests 
where the following sequence of cosmological epochs occurs: a. matter dominated phase (with or 
without usual matter), transition from decceleration to acceleration, accelerating epoch consistent 
with recent WMAP data b.ACDM cosmology without cosmological constant. As a rule, such mod- 
ified gravities are expressed implicitly (in terms of special functions) with late-time asymptotics 
of known type (for instance, the model with negative and positive powers of curvature). In the 
alternative approach, it is demonstrated that even simple versions of modified gravity may lead 
to the unification of matter dominated and accelerated phases at the price of the introduction of 
compensating dark energy. 

PACS numbers: ll.25.-w, 95.36.+X, 98.80.-k 



I. INTRODUCTION 

Modified gravity is extremely promising approach to 
dark energy. Within gravitational alternative for dark en- 
ergy (for review, see Q), the cosmic speed-up is explained 
by the universe expansion where some sub-dominant 
terms (like l/i?y,l3 or /ni?Q which may be caused by 
string/M-theory|3) may become essential at small cur- 
vature. It also explains naturally the unification of ear- 
lier and later cosmological epochs as the manifestation of 
different role of gravitational terms relevant at small and 
large curvature as it happens in the model with negative 
and positive powers of curvature @. Moreover, modified 
gravity may serve as dark matter. 

Special attention is paid to f(R) modified grav- 
ity which may be constrained from cosmologi- 
cal/astrophysical observational data and solar 
system tests |a la, S LLil- Recently, very interesting 
attempt to constrain such model (with positive and 
negative powers of curvature) from fifth force/BBN 
considerations has been made in 11] where it was shown 
that it is not easy to fulfil the known constraints and 
to describe the sequence of known cosmological epochs 
within the simple theory with positive and negative 
powers of curvature 6J. The cosmological dynamics of 
1/R and other related f(R) theories leading to late-time 
acceleration has been studied in refs.0,0,0,JEE3 while 
SdS black holes solutions were discussed in [la ]. 

In the situation when General Relativity can not nat- 
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urally describe the dark energy epoch of the univere the 
search of alternative, modified gravity which is consistent 
with Solar System tests/observational data is of primary 
interest. Of course, it is too strong (at least, at first 
step) to request that such theory should correctly repro- 
duce all known sequence of cosmological epochs (includ- 
ing inflationary universe where quantum effects may be 
essential). Nevertheless, it is reasonable to constrain such 
a theory (if this is really alternative theory for General 
Relativity) by the condition that it reproduces well estab- 
lished sequence of classical cosmological phases (matter 
dominated phase, transition from deceleration to accel- 
eration and current universe speed-up) being consistent 
with Solar System tests. In the present paper we con- 
struct several examples of such modified gravity where 
f(R) is presented implicitly, in terms of special functions. 

The important remark is in order. It is well-known 
fact that arbitrary f(R) gravity may be presented in 
mathematically-equivalent form as minimal scalar-tensor 
theory (Einstein gravity with scalar self-interacting field). 
Even more, it was shown HlHil that f(R) gravity may 
be formulated in mathematically-equivalent form as Ein- 
stein gravity with ideal field having inhomogeneous equa- 
tion of state. Then, it may look that it is not necessary to 
study modified gravity and it is enough to limit the con- 
sideration only by Einstein gravity with scalars and (or) 
ideal fluid. However, the situation is more complicated. 
For instance, modified gravity of specific form which de- 
scribes acceptable accelerating universe (with realistic 
effective equation of state) is not physically equivalent 
to scalar-tensor theory 0, 0]. Hence, the equivalent 
scalar-tensor gravity may not lead to accelerating FRW 
universe (or, it may lead but with significally different 
effective equation of state). The corresponding examples 
were given in [Til Il5|. Moreover, specific form of modi- 
fied gravity leads to specific form of scalar potential. As 
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a result, such specific modified gravity may comply with 
Solar System tests (acceptable Newton law, etc) while 
the corresponding scalar-tensor theory may not comply 
with it and vice-versa. Hence, one should consider all 
three classes of theories: modified gravity, scalar-tensor 
gravity and Einstein gravity with ideal fluid for descrip- 
tion of dark energy and early universe. One should fit 
all these three classes of theories with astrophysical and 
cosmological constraints (the corresponding cosm olog ical 
parameters are defined from different bounds in |l6jp in 
order to find finally what is the realistic gravitational 
theory compatible with observational data. In principle, 
In this work we are not so ambitious to comply with all 
bounds. We intend to present the first realistic exam- 
ple of modified gravity which is compatible with Solar 
System tests, which is cosmologically viable (sequence of 
matter dominated phase, transition from decceleration 
to acceleration, and acceleration phases) and which may 
lead even to ACDM cosmology. 



The paper is organized as follows. In the next section 
we present general formulation to reconstruct the mod- 
ified f(R) gravity for any FRW given cosmology (using 
the auxiliary scalar field). This formulation is applied 
to work out several models. The explicit example of the 
model where matter-dominated phase may be realized 
by pure f(R) gravity (no matter) with subsequent tran- 
sition to acceleration phase is presented. It is constructed 
the model where the function f(R) is expressed in terms 
of Gauss hypergeometric functions and where the stan- 
dard ACDM cosmology is reproduced. For the specific 
version of the modified f(R) gravity with matter it is 
shown that not only matter dominated phase with sub- 
sequent transition to acceleration occurs, but the accel- 
eration epoch complies with three years WMAP data. In 
order to ensure that transition from decceleration to ac- 
celeration indeed occurs, the (in)stability analysis of the 
cosmological solutions is fulfilled. We also demonstrate 
that corrections to Newton law for suggested versions of 
modified gravity are negligible. Third section is devoted 
to consideration of the modified gravity model [|| with 
compensating dark energy (ideal fluid). It is shown that 
the role of such compensating dark energy may be to 
ensure the transition from matter dominated to acceler- 
ation phase while during the current speed-up such com- 
pensating dark energy quickly disappears. Some outlook 
is given in the Discussion section. In the Appendix it 
is shown that our formulation is just equivalent to stan- 
dard metric formulation of f(R) gravity (without extra 
scalar). 



II. RECONSTRUCTION OF MODIFIED 
GRAVITY WHICH DESCRIBES MATTER 
DOMINATED AND ACCELERATED PHASES 

A. General formulation 

In the present section we develop the general formu- 
lation of the reconstruction scheme for modified gravity 
with f(R) action. It is shown how any cosmology may 
define the implicit form of the function /. The starting 
action of modified gravity is: 



S = 



-9f(R) ■ 



(1) 



First we consider the proper Hubble rate H, which de- 
scribes the evolution of the universe, with radiation domi- 
nance, matter dominance, and accelerating expansion. It 
turns out that one can find /(P)-theory realizing such a 
cosmology (with or without matter) . The construction is 
not explicit and it is necessary to solve the second order 
differential equation and algebraic equation. It shows, 
however, that, at least, in principle, we could obtain any 
cosmology by properly reconstructing a function f(R) on 
theoretical level. 

The equivalent form of above action is 



S = 



matter } • (2) 



Here P and Q are proper functions of the scalar field 
<t> and /Cmatter is the matter Lagrangian density. Since 
the scalar field does not have a kinetic term, it may be 



(con 

representation of f(R) gravity [lj|). In fact, by the vari- 

(3) 



regarded as an auxiliary field (compare with ideal fluid 
representation of /(it 
ation of 4>, it follows 



O = P'(#R + Q'(0) , 
which may be solved with respect to (f>: 

4> = 4>{R) ■ (4) 
By substituting (@J into (J2J, one obtains /(P)-gravity: 

S = J d 4 x^{f(R)+£ mllttcl .} , 

f(R) = P(^R))R + QWR)) . (5) 

By the variation of the action with respect to the 
metric g^, we obtain 



o = -^g^{P(4>)R + Q(m-R^P(<P) 
+V M V„P(0) - 5m „V 2 P(0) + 1t„„ . 



(6) 



The equations corresponding to standard spatially-flat 
FRW universe are 

r dP(<f>(t)) 



6H- 



-6P 2 P(0) - Q{cj>) 
Uh + 6P 2 ) P{<f>) + Q{cj>) 



dt 



P 



+2 



d 2 pm)) 

dt 



AH 



dPW)) 
dt 



P ■ 



(7) 



(8) 
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By combining © and Q and deleting Q(tf>), we find the 
following equation 



= 2 



d 2 Pm)) ^dPW)) 



dt 2 



-2H- 



■4HP(<f>) +p + p . (9) 



As one can redefine the scalar field <fi properly, we may 
choose 



t 



(10) 



It is assumed that p and p are the sum from the contri- 
bution of the matters with a constant equation of state 
parameters Wi- Especially, when it is assumed a combi- 
nation of the radiation and dust, one gets the standard 
expression 

p = p rQ a + p da a , p = — a , (11) 

with constants p r Q and p^o- If the scale factor a is given 
by a proper function g(t) as 



a = aoe 



s(*) 



(12) 



with a constant do, Eq.JHJ reduces to the second rank 
differential equation (see also 0): 

= 9 ^ ' d6 9 W> 

+ J2(l+w i )p i0 a^ 1+Wi) e- 3 ( l+ ^M . ( i3) 



In principle, by solving (|13|) we find the form of P{4>)- 
Using (JJJ (or equivalently JSJ), we also find the form of 
Q(4>) as 

Q{<P) = -6( 3 '(0)) 2 P(0)-6 ff '(^)^^ 

+ Y,P l oa Q 3{1+m) e^ 1+w ^ . (14) 

i 

Hence, in principle, any cosmology expressed as Ijl2(l can 
be realized by some specific /(-R)-gravity. 



B. Exactly solvable example I: unification of 
matter dominated and accelerated phases 

As an example, we consider the case 

g'{4>) = go + 9 i , (15) 

without matter p = p = for simplicity. Eq. (|13|l reduces 



n d2p ( ^9i\dP 2 9l 



whose solutions are given by the Kummer functions (hy- 
pergeometric function of confluent type) as 

P = z a F K (a, T ,z) , z 1 -TF K (a- 7 +l,2- 7 ;z). (17) 

Here 



z = go(p , a 



l + gi± y/gj + 2 5 i + 9 



7= 1± 



y/g( + 2 9l +9 



and the Kummer function is defined by 

"2, a (a + 1) • • • (a + n - 1) z n 



F K (a,j;z) 



7(7 + 1) ■ • • (7 + n — 1) n\ 



(18) 



(19) 



(20) 



Ea. (ll5|) tells that the Hubble rate H is given by 
H = g + T . 

When t is small, as H ~ gi/t, the universe behaves as 
the one filled with a perfect fluid with the EOS param- 
eter w = — 1 + 2/3<?i. On the other hand when t is 
large, H approaches to constant H — > g and the uni- 
verse looks as deSitter space. This shows the possibility 
of the transition from matter dominated phase to the ac- 
celerating phase (compare with 15]). We should note 
that in this case, there is no matter and /(i?)-terms con- 
tribution plays the role of the matter instead of the real 
matter. We will investigate later (in the next subsection) 
the example that there is a real matter. Similarly, one 
can construct modified gravity action describing other 
epochs bearing in mind that form of the modified gravity 
action is different at different epochs (for instance, in in- 
flationary epoch it is different from the form at late-time 
universe). 

We now investigate the asymptotic forms of f(R) in 
(J5J corresponding to l|15|) . When <f> and therefore t are 
small, we find 

P ~ P Q 4> a , Q ~ -6P <7i (gi + a) 4> a - 2 . 

Here Pq is a constant. Using 1(3}, it follows 

,2 6.9i (5i + a)(a- 2) 



(21) 



aR 

which gives 

f{R) „ _J^L ( 6gi (gi + a) (a-2) \ a/2 Rl _ 



(22) 



(23) 



On the other hand, when <fi and therefore t are positive 
and large, one gets 

P ~ P ^-^ ( 1 + (1 ' Q)(7 ' a) ) 



12^P o 2Q " 7 e 900 

-9 + 12a - 57 - 207 + 2a 2 



x 1 + 



25o <> 



9a 



,9o 



^-1 



V 12a 



(24) 
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Here Pq is a constant. Then we find 
1 

.9o 

R 



f(R) 



12.9o 2 



- 1 



9 „ 7 
h 9a 7 

2 2 ' 

-2Q+7+1 



2a — 7 



x exp 



9a 



R 



(25) 



This shows the principal possibility of unification of 
matter-dominated phase (even without matter!), transi- 
tion to acceleration and late-time speed up of the universe 
for specific, implicitly given model of f(R) gravity. 



C. Exactly solvable example II: model reproducing 
ACDM-type cosmology 

Let us investigate if ACDM-type cosmology could be 
reproduced by /(i?)-gravity in the present formulation. 

In the Einstein gravity, when there is a matter with the 
EOS parameter w and cosmological constant, the FRW 
equation has the following form: 



-H = p a 



-3(1+™) 



K 2 l 2 



(26) 



Here / is the length parameter coming from the cosmo- 
logical constant. The solution of l|26|) is given by 



a = doe 
9(t) 



g(t) 

2 Inf a sinhf 3(1 , + W) (t ~t ) j ).27i 



3(1 + to) V V 21 
Here to is a constant of the integration and 



1 



<* = -KfVpoa-o 



3(l+tu) 



(28) 



It is possible to reconstruct /(i?)-gravity reproducing 
1)2 7p. When the matter contribution is neglected , Ea. i|13|) 
has the following form: 

_^ sinh -(si±^ (( _ t0) )p M(29) 

By changing the variable from to z as follows, 
-2 f3(l + w) 



— — sinh 



2/ 



(t - 1 ) 



(30) 



Ea. (|29|l can be rewritten in the form of Gauss's hyperge- 
ometric differential equation: 



7 = 4 + 



d 2 P 
1 



= z(l-z)— + j- [a + /3 + l)z —-a(3P 



a/3 



3(1 + to) 
1 

"3(l + u;) 



5+^+1=6+ 



dP 

dz 
1 



3(1 + to) ' 



(31) 



whose solution is given by Gauss's hypergeometric func- 
tion: 



P 



P F{aJ,f,z) 



= Po 



T(a + n)T((3 + n) z r ' 



n=0 



T(7 + n) 



(32) 



Here T is the F-function. There is one more linearly 
independent solution like (1 — z)^ -0- ^^— a, J— (3, 7; z) 
but we drop it, for simplicity. Using 1)140. one finds the 
form of Q ((j)): 



P 

3(1 + w)z(l - z)Pg 

Z 2 (13 + 12w) 



F(a+l,/9 + l,7+l;«j|33) 



^From (|30|l . it follows z — * when t = <j) — * +00. Then 
in the limit, one arrives at 



P(<t>)R + Q(<P) 



6Pn 



Identifying 



Po = 



2k 2 



A = 



r- 



(34) 



(35) 



the Einstein theory with cosmological constant A can be 
reproduced. The action is not singular even in the limit 
of t — ► 00. Note that slightly different approach to con- 
struct ACDM cosmology from f(R) gravity is developed 
in ref.[l7|. 

Therefore even without cosmological constant nor cold 
dark matter, the cosmology of ACDM model could be re- 
produced by /(i?)-gravity. It was shown in ref.01 that 
some versions of modified gravity contain Big Rip singu- 
larities (for their classification, see Hjj). Hence, the 
above model without future singularity and with typical 
ACDM behaviour looks quite realistic. 



D. Models of f(R) gravity with transition of matter 
dominated phase to the acceleration phase 

Let us consider more realistic examples where the total 
action contains also usual matter. The starting form of 
g((j)) is 



g(</>) = h(</>)]n 



(36) 



with a constant cf>Q. It is assumed that h(<p) is a slowly 
changing function of cf>. We use adiabatic approximation 
and neglect the derivatives of h(<f>) (h'((f)) ~ h"(<f>) ~ 0). 
Ea. (|13|) has the following form: 



= 



d 2 P(<p) h{4>) dP (</))) 2h{<f>) 
d(j> 2 4> dt 2_ 

+E^ a o 3(1+u,,) ^ 3(1+tOi)s(0) 



P(4>) 



(37) 
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The solution for P(4>) is found to be 

P{(f>) = p+0 n+W) +p_0"- (<W 

+ J2^)^ 3{1+W ' M4,)+2 ■ (38) 



Here p± are arbitrary constants and 



ft(0) - 1 ± y/fogg + 6/i(0) + 1 
«±W = 2 ' 

/,x _ f/T . \ -3(1+™*) ,3ll+w)h(4>)\ 

PiW) = - {{l + w)pi a / 
x {6(1 + w){A + 3w)h(4>) 2 

-2(13 + 9w)h(4>) +4} _1 • (39) 

Especially for the radiation and dust, one has 



A J, 4h (<P) 



Pd{4>) = — 



3fl4 (40ft(</>) 2 - 32/i(<£) + 4) 



(40) 



a 3 (24ft(0) 2 - 26ft(0) + 4) ' 
We also find the form of Q(4>) as 

Q(4>) = -6H4>)P+ (h(<P)+n + (0))r +i,p) - 2 
-6ft(0)p_ (/i(0) + n_($) «/,"-^)- 2 

+ £){-6/i(0) (-( 2 + MK<t>) + 2) PiW 
+p. iO a o - 3(1+w) ^ (1+w) ' l(0) } 0- 3(1+w) ' i( « (41) 
Eq.® tells that 



H 



M*) 



and 



1? 



6 (— + 2ft(f) 2 



(42) 



(43) 



Let assume lim^o h(<f>) = hi and lim^oo h{4>) = hf. 
Then if < hi < 1, the early universe is in decceleration 
phase and if hf > 1, the late universe is in acceleration 
phase. We may consider the case h(4>) ~ ft-m is almost 
constant when <j> ~ £ m (0 -C t m <C +oo). If /ii, ft,/ > 1 
and < ft m < 1, the early universe is also accelerating, 
which could be inflation. After that the universe becomes 
deccelerating, which corresponds to matter-dominated 
phase with h(<f>) ~ 2/3 there. Furthermore, after that, 
the universe could be in the acceleration phase. 
The simplest example is 



K4>) 



hi + hfqtp 2 



(44) 



with constants hi, hf, and q. When ip — > 0, h(<p) — » hi 
and when (f> — > oo, h(cf>) — > hf. If q is small enough, 



h((f>) can be a slowly varying function of <f>. By using the 
expression of (|43|) . we find 



2 = $o(i?), *±(iJ), 



^ — 1/3 i 1/3 



a± 



A) 



2 (^2i? + 6ft/g- 12ftjg 

2i? + 6/i/g - 12ft 2 q) (ii + 6ft t g + 6ft/g - 4hih f q) 



+6ft, - 12ft 2 , 



3<?ii 



2R + 6h f q - Ylh f q 

3q 2 R 2 ' 
R + 6hiq + 6hfq — Ahjhfq 
q^R 



(45) 



There are three branches $o and <E>± in 1)45(1. Eas.(|43|) 
and 144fl show that when the curvature is small ((f) = t is 
large), we find R ~ 6 (—ft-/ + 2ft/) /0 2 and when the cur- 
vature is large = t is small), ii ~ 6 (—hi + 2hi) /<p 2 . 
This asymptotic behaviour indicates that we should 
choose $ m l|45fl . Then explicit form of f(R) could 
be given by using the expressions of P((f>) (|38(l and Q(4>) 
JUJ as 



f(R) = P (y$a(R)) R + Q (V$o(#) 



(46) 



One may check the asymptotic behavior of f(R) in 
(|46|l . For simplicity, it is considered the case that the 
matter is only dust (w — 0) and that p_ = 0. Then we 
find 



P(<f>) =p +< t> n+W +Pd 
One may always get 

n+ - (-3ft + 2) > 



3h(4>)+2 



(47) 



(48) 



in Here n+ = (h((f)) - 1 ± ^ft(0) 2 + 6ft(</>) + lj /2 

is defined in (|39|) . Then when <f> is large, the first term 
in l|47|) dominates and when 4> is small, the last term 
dominates. When <f> is large, curvature is small and 
dp 2 ~ 6 (—ft/ + 2ft/) /ii and ft(</>) — * ft(oo) = hf. Hence, 
Eq.||lZJ| shows that 



P(<t>)~P+ 



6(-ft/ + 2ft/)\(' l/ " 1+ ^ +6,! ' /+l)/4 



R 



and therefore 

/(ii) - R -(H4>)-5+^fh)+(,h } +i)/A 



(49) 
(50) 



6 



Especially when h ^> 1, we find 

f(R) ~ R- h f /2 . 



(51) 



Therefore there appears the negative power of R. As 
H ^ hf /t 1 if /if > 1, the universe is in acceleration phase. 

On the other hand, when curvature is large, we find 
2 - 6 (-hi + 2hi) /R and h(4>) -> h(0) = /i,. Then El 
shows 



P(4) ~ PdW 



-3hi+2 



(52) 



If the universe era corresponds to matter dominated 
phase (hi = 2/3), P(4>) becomes a constant and therefore 



f(R) ~ i? 



(53) 



which reproduces the Einstein gravity. 

Thus, in the above model, matter domi- 
nated phase evolves into acceleration phase 
and f(R) behaves as f(R) ~ R initially while 
f(R) ~ fl -("(«-5+V i J+ 6i '+ 1 )/ 4 at late time. 

Three years WMAP data are recently analyzed in 
Ref.[23|, which shows that the combined analysis of 
WMAP with supernova Legacy survey (SNLS) constrains 
the dark energy equation of state Wde pushing it to- 
wards the cosmological constant. The marginalized best 
fit values of the equation of state parameter at 68% con- 
fidance level are given by —1.14 < i«de < —0.93. In case 
of a prior that universe is flat, the combined data gives 
-1.06 < w D e < -0.90. 

In our model, we can identify 



WDE 



2 

3h~ f ' 



3 (1 + w D e) 



(54) 



(55) 



which tells h f should be large if hf is positive. For exam- 
ple, if wde = —0.93, hf ~ 9.51 • • • and if wde = —0.90, 
hf ~ 6.67- • •. Thus, we presented the example of f(R) 
gravity which describes the matter dominated stage, 
transition from decceleration to acceleration and acceler- 
ation epoch which is consistent with three years WMAP. 



E. (In)stability of the cosmological solutions 

Let us investigate the stability of the obtained solu- 
tions. We assume 



a = a e gW 



(56) 



which corresponds to l|12|) and P(<j>) should be given by 
a solution of l|13|) (and Q(<p) should be given by 11411 ). 



Under the above assumptions, we consider the perturba- 
tions in Q) and ©. By deleting Q((j>) in Q and ©, one 
obtains 



= 2 



d 2 P(0) „ ,,^dP(4>)) 



dt 2 



2g'(4>)- 



dt 



4<7"(</>)P(0) 



■J2(l + Wi )p 



d 2 P((/)) (d<t> 



d(j) 2 



dt 



dP(4>) d 2 (j) 
dd> ~dl 2 



y yv) dd> \ dt 



+ 2(l + u,) M a - I(,+ -')e- 3 ('+-W*> . (57) 



Then combining ljl3|l with Q56JI it follows 



= 2 



d 2 P(<j>) gf{(f)) dFm 




2</'(</>)P(</>) 



By defining 5 as 



dt 



(58) 



(59) 



we consider the perturbation from the solution l|10|) . Us- 
ing (J5SJ, one gets 



(60) 



dP(4>) 



0£l +2 g>(t)P(ct>) 



Then when lj > (to < 0), the perturbation becomes 
small (large) and the system is stable (unstable). 

As an example, the case <|36fl may be considered. It 
gives (|3"5j) with J2U). Especially when p± = p i0 = 
except pdi, we find 



U) = 



2(l2/i 2 - 13/i + 2) 
(2 - h)t 



(61) 



Here the derivatives of h(cf>) like h'((f>) are neglected again. 
Then u) goes to infinity when h = 2 and h — ► ±oo and lo 
vanishics when 



(62) 



h = h±= 13 ^7 63 = 0.87238 • • ■ , 0.2109478 • • ■ 

24 
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Therefore uj > and the system is stable when ft < h— 
or h + < ft < 2. Hence, for the case that the universe 
starts from decceleration phase with h = ho < 1, if ho > 
ft + , there is a stable solution where universe develops 
to the acceleration phase h — *■ 2 > 1. Even if we started 
with h(4>) = 2/3, which corresponds to matter dominated 
phase a ~ t 2 ' 3 , the solution is unstable since ft_ < 2/3 < 
ft + , the perturbed solution could develop into the stable 
solution with h > h + and therefore there could be a 
transition into the acceleration phase. If h goes to 2 
from the region with h < 2, since uj — > +oo, the solution 
becomes extremely stable. Hence, h may pass through 
the point h — 2 and h could become larger than 2, where 
the effective EOS parameter w — 1 — 2/3ft = 7/9. 
Note that when p± = pio = (i ^ d : dust), one gets 



In the expression l|67[l , uj vanishes when h — h± (|62|l and 
h = h±, defined by 

h± = 4± x ^ = °- 6449 ' • • : 0.15505 • • • . (69) 
Under the assumption JfJSJ), uj diverges at 



R 



f(R) ~ R 3h ' 2 



(63) 



Therefore in the matter dominated phase h ~ 2/3, the 
action behaves as the Hilbert-Einstein action. 

In more general case that there is only one kind of 
matter with w and p± = 0, we find 



R 



3ft(</>) {12 (1 + w) h((f>) 2 - 2 (7 + 9w) h(cf>) + 4} 



{3(l + w)h((b)+2} <j> 2 
f(R) ~ rH^WHR)) _ 

Note that to in lj60|l is given by 

_ 3(1 + w)(4 + 3w)ft(<?(>) 2 - (13 + 9w)/t(</>) + 2 
{-(l + 3w)ft(0) + 2}/ 



(64) 



(65) 



3a o a o 
x (12ft 2 - 13/i + 2) 



3a o 

3fc 



^10ft 2 -8/i + l) (ft -2) I 



(67) 



As clear from lfTT|) . iJSSJ, and p r o3aQ 4 and /Odo^o 3 

correspond to the energy density for radiation and dust 
(usual matter plus cold dark matter), respectively, when 
t = <p = 4>q. Let us choose t = <f>o corresponding to the 
present universe. It may be assumed 



3<4 



_ f±p r o<t>f<t> h \ ( Pd o4>l h 



\ n 3 
\ "0 



< 1 



(68) 



24 
25 f 



S± 



(28+17V6) (-6±V6)(l6±x/6) 
+ 25000 ■ 1 ' 



Here <$± > 0. Near the singularities, if uj is negative, there 
could be very large instability, which should be avoided. 
In case p± ^ 0, the singularity could be avoided. The 
points where uj vanishes could remain even if p± =/= 
but the instability becomes finite and there can be a 
solution which describes the transition from the matter 
dominated phase to the acceleration phase. 

We now consider the case that the contribution from 
the matter could be neglected in l|38|) and therefore we 
could assume that p^ = 0. Furtheremore when one of 
p± vanishes, one finds 



R 



f(R) ~R N ±, N± = 1 - '-f . (71) 



The behavior of N± is as following 



Next we consider the case that pdo 7^ and p r o 7^ 
but p± = and p^ — except i 7^ d, r. Then it follows 

ui = -i {4p r (I0h 2 -8ft+ 1) + 2p d (I2h 2 - 13/i + 2) <j) h ) 

x {2p r (h-l)+p d {h-2)cb h y 1 , (66) 
or by using (|39"|) . 

2/W , P ^\, mh2 _ 8h+l) 



when ft -> , N + -> 1 and iV_ 



ft -> +00 , 2V+ -»• -- and N- 
ft = 1 , 7V± -> 1 T 



V2 



ft = f , ^4 



1 and iV_ - I . (72) 



Thus, in case p+ = but p_ 7^ 0, when ft, — ► +00, the 
higher derivative inflationary model f(R) ~ i? 2 appears. 
We should also note w in (|60|) has the following form: 



uj = UJ± 



2 (ft(0) - 1) ± Ayjh{4>) 2 + 6ft(</>) + 1 
5ft(</>) - 1 ± ^h{<j)) 2 +6ft(0) + 1 



(73) 



Now ft± may be defined as 



ft± 



13 + 4V10 



8.5497- ■ -,0.11690- ■ 



(74) 



Then uj + < when ft > ft + or ft < 0, and lj + > when 
< ft < ft + . We also find that when cj_ > 0, ft > 2/3 
or ft < ft_, and lj_ < when ft_ < ft < 2/3. Hence, the 
model where p+ = but p_ 7^ is stable when ft > |. 
Thus, one can make a stable model where matter dom- 
inated phase ft ~ 3/2 evolves to the acceleration phase 
ft > 1. 
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We should note that p- = and pio = but p+ ^ 0, 
if we put 



h 



10 

y ' 



we obtain N + = — 1 or f(R) ~ 1/i?,. Let us put 



10 



It follows 



/i(0) = — + Sh , \Sh\ < 1 



1 R 

N+ = -1 - — <5/i 
17 



(75) 



(76) 



(77) 



When ft ~ 10/3, from ©, the curvature is given as 



R 



6h(h + n + )(n + - 2) 



n +4> 2 



220 

302 



(78) 



Hence, with the choice 



Sh 



17 



220 2 



220 3 /? 



17 



R 2 pR 3 



18 In 4 V 2 «V 



(79) 



with a constant /z, which has a dimension of mass, one 
arrives at 



fiR) R + 2^ 



(3R 2 



(80) 



which reproduces the action proposed in 6|. (Note that 
such class of actions does not describe sequence of matter 
dominated/acceleration phase H^)- Thus, using stability 
analysis we demonstrated that indeed the matter dom- 
inated phase may transit to the acceleration phase for 
some implicit model of f(R) gravity found in the previ- 
ous subsection. Moreover, it is shown that the model of 
ref. with positive and negative powers of the curvature 
is just asymptotic form of such consistent f(R) theory (at 
some specific values of parameters) at acceleration epoch. 
The complete, implicit version of f(R) theory found in 
previous subsection decribes the sequence of matter dom- 
inated phase, transition from deceleration to acceleration 
and then acceleration epoch of the universe. 



F. Corrections to Newton law 

In the present subsection we will discuss the contribu- 
tions to Newton law in the modified gravity under consid- 
eration. Note that there is big number of papers devoted 
to the study of newtonian regime in modified f(R) grav- 
ity 0,0- However, these papers are mainly devoted to 
the study of newtonian regime for 1 / R models. We now 
check when the corrections to the Newton law are not 



essential in /(i?)-gravity under consideration. For this 
purpose, we put a point source at r = 0: 



m xr \ 



a(tf 
Transforming 



(81) 



(82) 



one finds the (t, £)-component of @ has the following 
form: 

= Go + G\ , 
Go = ~Sgu {P{<t>) (pH + 12tf 2 ) + Q(0)} 

~P{<t>)8R-P{ff>)8Rtt , 
Gi = i{p'(0)(6ij + 12ff 2 ) +Q'(0)}ty 
+3 (h + H 2 ) P' 



--P'{<t>)V t 8gu + 3HP'(cf>)5g tt + P"^)Sg tt 



-HP'(<f>)g^S 9ij + \p'{<t>)Vt (g^Sg^) 
1 -S(r) . 



(83) 



2a{tf 

Here the gauge condition is chosen 

VSg^ = . (84) 
On the other hand, Eq.© gives 

= {P"{<j>) + Q"(4>)R) 5<j) + Q'(4>)SR . (85) 
We now consider the region for r = \r\ as 



or 



1 1 

- « r < — , 
m H 



m > — » H 
or 



(86) 



(87) 



We should also note that 



Then if 



P'{4>) < HP{(f>) , P"{(f>) < H 2 P{(f) , (89) 

Gi (|%3f could be neglected if compared with Go- Then 
since <j> can be regarded as a constant and therefore 8<fi = 
as long as is satisfied, Ea. (|53l reduces to that in 
the Einstein gravity by identifying P(<fi) with 1/k 2 and 
Q(4>) — A 2 /k 2 since P(4>) and Q{<f>) are slowly varying 
(almost constant) functions. From l|85|) . 5<fi — implies 
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SR = 0. Then Ea. (|83(l gives (t, i)-component of the usual 
perturbation in the Einstein equation: 

= G 

= -^9tt{Ro+A 2 }-^(^SR-8R t 

R a = 6H+ 12H 2 . (90) 

Here clq is the scale factor in the present universe ao = 
a(i), which may be chosen to be unity ao = 1. For the 
almost flat universe as our current one, we can neglect 
the first term in (|§U)l : R + A 2 ~ and we obtain 



= -4(V 



SR t 



md(r) 



(91) 



which further reduces, since the universe is almost flat 
and the source and the universe are static in a region 
where we are investigating the Newton law, to 



0=2^(A<5 5ft 



(92) 



Here A is usual Laplacian. Therefore if the conditions in 
(|88|l are satisfied, the correction to the Newton law could 
be neglected. In case of (|42|) . if we choose p+ and p- 
to satisfy (|S5|) in the present universe, the correction to 
the Newton law could be small. (It is interesting to note 
that standard Newton law is valid also in an arbitrary 
F(G) gravity [23| where G is Gauss-Bonnett invariant). 
One may consider even simpler situation: admitting that 
Newton law is satisfied only in current universe. Then, 
the form of Newton law should be fixed only at present 
universe (in other words, the initial value of f(R) should 
be restricted). We should also note that if P(cj>) changes 
its sign, as we identified P(<fr) with 1/k 2 , there could 
appear the anti-gravitation effect. For instance, the form 
of modified gravity may be changed in the future, at the 
end of acceleration epoch, driving the Newton law to its 
opposite sign form. 



III. MODIFIED GRAVITY AND 
COMPENSATING DARK ENERGY 

In the present section we will present another approach 
to modified gravity. Specifically, we discuss the modi- 
fied gravity which successfully describes the acceleration 
epoch but may be not viable in matter dominated stage. 
In this case, it is demonstrated that one can introduce 
the compensating dark energy (some ideal fluid) which 
helps to realize matter dominated and decceleration- 
acceleration transition phases. The role of such compen- 
sating dark energy is negligible in the acceleration epoch. 

We now start with general /(i?)-gravity action: 



S = J d i x{f(R) + C 
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(93) 



In the FRW metric with flat spatial dimensions one gets 



P = f(R) -6[H + H 2 - H- f(R) , 
p = -f(R) - 2 (-H-3H 2 + ^ + 2H^ ) f(R) , 



dt 



R = 6H + 12H 2 



(94) 



If the Hubble rate is given (say, by observational data) 
as a function of t: H = H(t), by substituting such ex- 
pression into H94|) , we find the ^-dependence of p and p as 
p = p(t) and p = p(t). If one can solve the first equation 
with respect to t as t = t(p), by substituting it into the 
second equation, an equation of state (EOS) follows: 



P = p(t(p)) 



(95) 



Of course, p and p could be a sum with the contribution 
of several kinds of fluids with simple EOS. 

We now concentrate on the case that f(R) is given 

b y 



(3R 2 



(96) 



Furthermore, we write Hit) as Eq. I|42(l and assume h{t) 
is slowly varying function of t and neglect the derivatives 
of h(t) with respect to t. Then one gets (|43|l . 

First we consider the case that the last term in H80JI 
dominates f(R) ~ PR 2 , which may correspond to the 
early (inflationary) epoch of the universe. It is not diffi- 
cult to find 



P 



36/3(~l + 2/t(<))/i(t) 2 
36/3(-l + 2/i(f))/j(3/i(i) + l) 



(97) 



If h goes to infinity, which corresponds to deSitter uni- 
verse, we find p ~ p ~ h 3 although, from (|43|) . R ~ h 4 . 
Therefore p, p <g; (3R 2 and contribution form the matter 
could be neglected. Then the inflation could be gener- 
ated only by the contribution from the higher curvature 
term. 

Second, we consider the case that the second term in 
(|SU|l dominates f(R) ~ A, which may correspond to 
the matter dominated epoch after the inflation. In this 
case p and p behave as 



P 
P 



I2h(t) + 6h(t) 2 
^H 2 ' 
4h(t) - 6h(t) 2 
^t 2 



In the matter dominated epoch, we expect h 
a ~ ti ) . Hence, one gets 



32 
3^ 



p~0 



2/3 



(99) 
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Therefore in the matter sector, dust with w = (p = 0) 
should dominate, as usually expected. 

Finally we consider the case that the first term in 1(93() 
dominates f(R) ~ —a/R n , which might describe the ac- 
celeration of the present universe. The behavior of of p 
and p is given by 



p ~ a {6 (n + 1) (2n + 1) h(t) + 6 (n - 2) h(t) 2 } 
x {-6ft(£) + 12ft(£) 2 }~"~V™ , 

{-4n(n + l) (2n + 1) - 2 (8 ra 2 + 5n + 3) ft(t) 
-6(n-2)ft(i) 2 } 

x {-6h(t) + 12fc(f) 2 }~ n_1 t 2 ™ . (100) 
Thus, the effective EOS parameter wi is given by 



p 
p 



~ {-An (n + 1) (2n + 1) - 2 (8n 2 + 5n + 3) ft(£) 
-6 (n - 2) h(t) 2 } {6 (n + 1) (2n + 1) ft(£) 
+6(n-2)ft(i) 2 }~ 1 . (101) 

In order that the acceleration of the universe could occur, 
we find ft > 1. Let us now assume that h(t) — > ft/ when 
£ — * oo. Then one obtains 

Wt —> Wf 

= {-4n(n+ 1) (2n+ 1) - 2 (8n 2 + 5n + 3) ft/ 
-6 (ti - 2) ft 2 :-} {6 (n + 1) (2n + 1) ft/ 
+6(n-2)ft 2 } -1 , (102) 

and -ff(£) — > ft/A- Since the matter energy density p Wf 
with the EoS parameter Wf behaves as 



p Wf a 



-3(l+i"/) 



oc 



exp(-3{l + w f ) J dtH(t)\ 



the energy density is 

p Wf oc t -^ 1+w ^ h f . 
Comparing 1 104(1 with il 100(1 , we find 

2n = —3(1 + wf)hf , 



(103) 
(104) 

(105) 



which can be confirmed directly from (|102(l . 

^From the above consideration, we find p and p contain 
mainly contributions from dust with w = 0, pd{t), p(t) = 
and "dark energy" with w — wi in (|101(l . p;(t). 
In the expressions of p(t) and p(t) in l|94|) . there might 
be a remaining part: 

PR (t) = p{t)-p d {t)-pi{t) , PR {t) EE p{t)-pi{t) , (106) 

which may help the transition from the matter dominated 
epoch to the acceleration epoch. By deleting t in the 



expression of 1(106(1 . we obtain the EOS for the remaining 
part: 



Pb. = Pr(pr) 



(107) 



which may be called the compensating dark energy. More 
concretely, according to 1(99(1. one may have 



32 



- 3 r* dtiM 
j t t 



(108) 



and according to 1(100(1 . 

pi - a{6(n+l)(2n + l)ft/ + 6(n-2)ft}} 

x {-6ft/ + 12ft 2 if e - 3(1+ ^ dt ^fl09) 

In ((109(1 . we choose ii to be large enough. When t ~ t , 
p(i) ~ and when £ — > oo, p(£) ~ p;. Thus, only 
dominates after t = t\ but it becomes smaller in late 
times. Hence, the role of pn (which perhaps may be 
identified partially with dark matter) is only to connect 
the matter dominated epoch to the acceleration epoch. 

IV. DISCUSSION 

In summary, we developed the general formulation of 
modified f(R) gravity which may be reconstructed for 
any given FRW metric. The resulting action is given in 
the implicit form (usually, in terms of some special func- 
tions). Nevertheless, its early and late times asymptotics 
may be defined, it turns out to have quite simple form, 
for instance, as model with negative and positive pow- 
ers of curvature. Several examples predicted by realistic 
cosmology are constructed. These specific f(R) theo- 
ries describe the sequence of cosmological epochs: mat- 
ter dominated stage (if necessary even without matter), 
transition from deceleration to acceleration and cur- 
rent cosmic speed-up consistent with three years WMAP 
data. Moreover, the study of their newtonian regime 
indicates that such models are consistent with Solar Sys- 
tem tests. It is not difficult to extend such a formulation 
to include consistently also radiation dominated phase 
(perhaps, even inflation). Hence, modified f(R) grav- 
ity indeed represents the realistic alternative to General 
Relativity, being more consistent in dark epoch. It is 
also shown that some implicit version of f(R) gravity 
may describe ACDM cosmology without need to intro- 
duce the cosmological constant and without singularity 
near R = 0. 

In the alternative approach we also demonstrate that 
even simple models like the ones of ref.0, become 
cosmologically viable if compensating dark energy is in- 
troduced. It remains to study if such compensating dark 
energy or the version of f(R) gravity which mimics mat- 
ter dominated phase without matter may serve as dark 
matter of the universe. 

Definitely, more careful study of modified gravity 
and fitting the above models against the observational 
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data/various constraints [T|j should be done. First of all, 
one should study the linear perturbations in the matter- 
dominated epoch. Such a study made for scalar-Gauss- 
Bonnet gravity model of ref.|2J| in ref.[25| shows that 
there is no really strong change if compare with usual 
General Relativity. Hence, one may expect that the same 
will occur in the present model while careful study should 
be made, of course. It is also known [jj that Supernovae 
la constraints are easy to fit in the class of models un- 
der consideration (subject to assumption that they are 
treated as usual candles). Second, the CMBR peak loca- 
tions are weakly model-dependent. Nevertheless, it is im- 
portant to check constraints appearing from the CMBR 
shift parameter and baryon oscillations as well as nu- 
cleosynthesis bounds which restrict the amount of dark 
energy in the current universe. Third, Solar System con- 
straints (time variation of the effective gravitational con- 
stant, study of PPN parameters) should be considered 
in detail. The preliminary expectation is that this may 
be achieved due to the fact of freedom of the form of 
modified gravity (as well as first time derivative of the 
gravitational Lagrangian) at some specific (initial) time. 
This will be investigated in detail elsewhere. Having in 
mind that new, more precise observational data will be 
available soon, one may expect that the question: is mod- 
ified gravity suitable as dark energy will be answered in 
near future. 
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APPENDIX A 

Let us show for the explicit example of f(R) that our 
formulation with auxiliary scalar fixed as time is equiva- 
lent to usual metric formulation. The starting action of 
the modified gravity coupled with matter is: 



S= / d 4 x^{f a R a +£ maUer } 



(Al) 



FRW equation is given by 



1 



= fo< -§ (6H+12H 2 



-3a (H + H 2 ^ (qH+12H : 



(A2) 



An exact solution of (|A2|) is given by 

h , 2a 

a = a t , ho 



3(1 +w) 



ao 



Pa 



x {(1 - 2a) (1 - a) - (2 - a)h }]~^T^) . (A3) 

Instead of action IjAljl . one now starts with the action 
©• Eq.JXSJl shows 



H = — or g(t) = ho lnt . 



(A4) 



Hence 




?"(0)P(0) 



+ (1 + W ) poflo 3(1+tu) e- 3(1+t0)9(t) 
d 2 P{4>) 2h dP{4>)) 4^o 



dcj> 2 <j) d(j) <f 
+ (l + w)p Q a '4> 



P(<t>) 



(A5) 



Here we have used a relation h — (2/3) (a/(l + w)) in 
(|A3() . A solution of l|A5|l is given by 



Po 



Then we find 



P = P o 0- 2a+2 , 

[l + w)p a 



4{(l-a)(l-2a)-/i (2-a)} 



(A6) 



Q(4>) = -6(.g'(0)) 2 P(0)-6 g '(0)^M 

+p a- 3il+w) e-^ 1+w ^ 

Q = -6h Q (h a -2a + 2)P + p a- 3il+w) 
6(1 - a)h (2h - 1)P 



In the last line, we used the definition of Pq in (|A6|) 
Therefore it follows 



(A7) 



= P'{4>)R + Q'{4>) 
which gives 



= 2(1 - a)<i>- 2a+1 R - 2aQ (j>- 2a , (A8) 



2 _ aQ I2hl - 6h 



(1 - a)P R 



R 



(A9) 
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Then the action© has the following form 

S = J d A X^g{foR a +C ma tter} 

fo =- 



1 — a 



—ot+l f^—a+l pa 
^0 M) 



-3(1- 



(12/iq - 6h ) a poa 
~ 6{(l-a)(l-2a)-ft (2-a)} 

Comparing f' i|A10|) with l|A3(l . one gets 
fo = fa i 



which shows the action lAlj l is surely reproduced. On 
the other hand, Ea . (|A3|1 shows 



R= P 



(A10) 



(All) 



Comparing (|A12|) with (|A9|) , we find 



(A12) 



(A13) 
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